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We present the temperature and density dependence of symmetry energy for nuclear matter at finite temperature based on the 
approach of the thermodynamics with Skyrme energy density functional. We first classify the Skyrme interactions into 7 groups 
according to the range of neutron and proton effective mass in neutron matter limit(99.99 per cent neutron in the matter). We find 
that there is obvious correlation between the temperature dependence of the symmetry energy and the splitting of the neutron and 
proton effective mass. For some Skyrme interactions with m* > m* and strong splitting of the neutron and proton effective mass 
in asymmetric nuclear matter, a transition of the temperature dependence of symmetry energy from decreasing with temperature 
at low densities to increasing with temperature at high densities appears. For other Skyrme interactions, we do not observe such 
phenomenon. Our study show that the symmetry energy in hot asymmetric matter not only depends on symmetry potential part but 
also on the splitting of the neutron and proton effective mass to a certain extent. 

Keywords: nuclear symmetry energy, nuclear symmetry energy at finite temperature, Skyrme density functional, Splitting of 
neutron and proton effective mass 



The symmetry energy for nuclear matter iisymis ver y impor- 
tant for understanding not only the structure of nuclei far from 
8 stability line, the dynamics and many interesting phenom- 
ena in reactions of isospin asymmetric nuclei but also many 
critical issues in astrophysics! ill- Recent research on the sym- 
metry energy for cold nuclear matter has already obtained a 
great progress and some constraints on the density dependence 
of the symmetry energy at subnormal densities have been set 
|S S B S 01 • But the density dependence of the symmetry 
energy at high densities extracted from the ratio of 7t~/tt + and 
other observables in heavy ion collisions at intermediate high 
energies are extremely model dependent and divergent Si- 
Moreover, for heavy ion collisions at intermediate high ener- 
gies, the system should be heated and therefore in addition 
to the density dependence, the temperature dependence of the 
symmetry energy is also of fundamental importance for heavy 
ion collisions with nuclei away from /^-stability line as well as 
for astrophysics. While, how the symmetry energy depending 
on temperature is even unclear H9l 1 1 2L 1 1 Oil and only few informa- 
tion on the symmetry energy at very low density and very low 
temperature was extracted from experimental measurement in 
heavy ion collisions up to now Bl ill 

The symmetry energy for hot nuclear matter can be expressed 

as 



E sym (p, T) = E(p,T,S=l)- E(p, T, 6 = 0). 



(1) 



E(p, T, 6) is the energy per nucleon for nuclear matter with den- 



and p are the neutron, proton and nucleon density of the system, 
respectively. It gives an estimation of the energy cost to con- 
vert all protons in symmetric nuclear matter(SM) to neutrons at 
fixed temperature T and density p. 

In addition to the symmetry energy, the splitting of the neu- 
tron and proton effective mass (EFM) is also a hot topic in the 
study of the properties of asymmetry nuclear matter ll 1 3l 1 1 411 . 
The neutron and proton EFM splitting predicted by different ef- 
fective interactions are rather different^ [3 O ES E • 
Both the symmetry energy and the neutron and proton EFM 
splitting closely relate to the isospin dependence of the nuclear 
interactions. Then one would naturally ask whether the splitting 
of the neutron and proton EFM influences the temperature and 
density dependence of the symmetry energy? In this letter we 
will first study the symmetry energy at finite temperature and 
then investigate the correlation between the EFM splitting and 
the density and temperature dependence of symmetry energy. 
The approach of thermodynamics with Skyrme energy density 
functional is applied in this work. 

The Skyrme energy density functional reads 

H = £[ T/ ,(r) + r„(r)] 

+|f [(2 + x )p 2 - (2x + l)(p 2 p + P 2 „)] 
+±hp a [{2 + x 3 )p 2 - (2x 3 + l)(p£ +pl)l 
+|[?i(2 + xi) + f 2 (2 + x 2 )]rp 
+|[f 2 (2x 2 + 1) - h{2x\ + 1)](t pPp + r„p„), (2) 



sity p, temperature T, isospin asymmetry 6- 



_p,,-Pp 



where p„, p p , where p = p„ + p p , and r = t„ + t p . The p q and T q (q denotes 
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protons or neutrons) are computed by 
n q (p)^§~dp, 

o 



and 



/-*co 

= 2 I n q (p)^f^dp. 
Jo 



(3) 



(4) 



The occupation number distribution for species q, n q (p), obeys 
the Fermi-Dirac distribution function, which reads 



n q (p) = 



1 



1 + exp|j6(e 9 - fi q )] 



(5) 



sjjjjjr + I/ 9 is the single particle energy for species q. The 



EFM m* for species q is expressed as 



2m- 



(6) 



^ + \[h(2 + x { ) + t 2 (2 + x 2 )]p 
+\[h(2x2 + Y)-h(2x\ + l)]p q . 
The mean field U q for species q reads 

U q (r) = if [(2 + x )p - (1 + 2x )p q ] 
+ ±t 3 [(2 + x 3 )(2 + a)p a+1 
-(1 + 2x i )[2p a Pq + ap^Upl+pj)]] 
+g[fi(2 + xi) + t 2 (2 + x 2 ~)]T 
+ i[f 2 (2x 2 + l)-h(2 Xl + l)] Tq . (7) 

Introducing effective chemical potential p q = p q — U q , then 

= ; . L : — • (8) 



1 + exp[J3(p 2 /2m* - p' q )]' 

By solving equation © and (JHJ with EFM depending on density 
iteratively, for any pair of p n and p' p we can obtain the proton 
and neutron density p p and p„ at temperature T. Then, the en- 
ergy per nucleon in neutron matter(NM) and SM can be calcu- 
lated by using equation (0]i and ^ and finally symmetry energy 
at finite temperature can be calculated. 

Since the nucleon EFM is explicitly involved in the itera- 
tion process we first make a general survey of the density de- 
pendence of the neutron, proton EFM for various Skyrme in- 
teractions. 94 Skyrme interactions are tested and they can be 
divided into 7 groups according to their corresponding ranges 
of the neutron and proton EFMes varying with density. Ta- 
ble Q] presents the groups of Skyrme interactions and the cor- 
responding ranges of m*/m in SM and NM. The m is the bare 
nucleon mass. In this work, the NM is always referred to a 
nuclear matter limit with 99.99 percent neutrons in the mat- 
ter. For convenience, the EFM m* discussed in the following 
text always refer to the ratio m* jm (the scaled effective mass), 
which should not be confused with the EFM m* appeared in ex- 
pression (6) and (8). In SM, the neutron and proton EFM are 
equal and all Skyrme interactions are simply in the groups: (I) 
m* > 1 and it increases with density increasing; (II) m* = 1 and 
it is independent on density; and (III) m* < 1 and it decreases 
with density increasing. Most of Skyrme interactions belong 
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Figure 1: (color online) Density dependence of symmetry energy at T=0, 5, 
10, 20 MeV calculated with different Skyrme interactions. The ranges of cor- 
responding neutron and proton EFM and the ratio between the neutron EFMes 
in NM and SM at saturation density are also presented in each sub-figure. 



to the group (III). For asymmetric nuclear matter, the situation 
becomes more complex because of neutron-proton EFM split- 
ting. The splitting of the neutron and proton EFM in asymmet- 
ric matter is currently not known empiricaly[9]. Theoretical re- 
sults on the splitting of the neutron and proton EFM are highly 
controversial among different approaches and different effective 
interactions. One can see from Table I that only few Skyrme 
interactions belong to groups 1,11, and III and most of Skyrme 
interactions belong to groups IV to VII with splitting of neutron 
and proton EFM. Large part of them belong to groups IV and 
V (m* > m* p ), others belong to groups VI and VII (m* < m* p ). 

Now let us come to study the density and temperature de- 
pendence of the symmetry energy in hot nuclear matter. We 
show in Fig. Q] the symmetry energies in nuclear matter at 
temperatures T- 0, 5, 10, 20 MeV calculated with different 
Skyrme interactions, respectively, as typical examples. The 9 
subfigures in FigQ]is ordered according to the magnitude of R® n 
for the Skyrme interaction applied. Here we introduce the ra- 
tio R m = t-, where the subscripts and 1 indicate the isospin 
asymmetry 6 =0 (for SM) and 1 (for NM), respectively, to 
characterize the strongness of the splitting of the neutron and 
proton EFM for the Skyrme interaction applied. The quantity 
R® n presented in each subfigure is the value of R m at normal 
density. Obviously, the R m is proportional to the neutron and 
proton EFM splitting in asymmetric matter. The values of R m 
for Skyrme interactions in groups IV and V are always larger 
than 1 and those in groups VI and VII are smaller than 1 . As 
is known that the symmetry energies for cold nuclear matter 
calculated from different Skyrme interactions are largely di- 
vergent especially at high densities. We can also classify the 
Skyrme interactions into two groups A and B according to the 
trend of the density dependence of the symmetry energy in 
cold matter. For the group A, namely the Skyrme interactions 
like SLy7, SkT5, SkMP, and SkM, etc., nuclear symmetry en- 
ergy increases monotonously with density and for the group B, 
namely the Skyrme interactions like Skz3, SHI, SKXm, SkP, 
and v070, the nuclear symmetry energy first increases with den- 
sity until certain density then it bends down as density further 
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Table 1: The ranges of the neutron and proton EFMes, m* /m, in SM and NM calculated with various Skyrme interactions. 
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increases and finally becomes negative. The group A corre- 
sponds to the group I and the group B corresponds to group 
II and III in ref . 12111 . respectively. It is seen from Fig.l that 
there is no obvious correlation between the trend of the den- 
sity dependence of the symmetry energy in cold matter and the 
magnitude of the R® n of the corresponding Skyrme interactions. 
As the temperature dependence of the symmetry energy is con- 
cerned, one sees that the nuclear symmetry energy decreases 
with temperature increase in all densities for Skyrme interac- 
tions with small R° M , such as SLy7, Skz3,Skt5,SkMP,etc, in 
consistency with the results given in [9]. But for SkM, SKXm, 
SkP, and v070 with large R^ the symmetry energy decreases 
with temperature increasing at low density and increases with 
temperature when the density is higher than a certain density. 
We call this phenomenon the transition of temperature depen- 
dence of the symmetry energy(TrTDSE). From Table[T]one can 
find that Skyrme interactions such as SkM, SKXm, SkP, and 
v070 for which the TrTDSE phenomenon appears, all belong to 
groups IV and V satisfying m* > m* p . The density for the on- 
set of the TrTDSE depends on the magnitude of the R^, of the 
corresponding Skyrme interactions. The larger R 1 ^ is, the lower 
density is for the onset of the TrTDSE. For the case of SIII, 
we find that the symmetry energy at finite temperature is very 
close to that of cold matter up to p < 4po where our calculation 
stops. For Skyrme interactions belonging to the groups I-III and 
VI, VII (not satisfying m* n > m* ), the TrTDSE phenomenon can 
not happen. 

In order to explore the condition of TrTDSE in hot matter 
more quantitatively we study the difference between the sym- 
metry energy in a matter at temperature T and that in a cold 
matter, which reads as 

AE sym (p, T) = £ sym (p, T) - E sym (p, 0) 

= (£+*3+*4)(tT-t?)(1-£), (9) 

with b 3 = |[fi(2 + xi) + f 2 (2 + x 2 )] and b 4 = + 2x 2 ) - 

f i ( 1 + 2x\)]. It is easy to extend equation (9) to E sym (p, AT) 



with AT = Ti — T\ , The R T in (9) is the ratio between the ki- 
netic energy density increasing from T = to T = T in NM 

T T -T° 

and that in SM, R T = ' T _ ' , where the subscripts and 1 in- 
dicate the isospin asymmetry 6 =0 and 1, respectively. The 
first factor in the expression (|9]l depends on the Skyrme interac- 
tion parameters which is always positive up to p = 2po for the 
94 Skyrme interactions involved in this work( in fact for some 
Skyrme interactions this density can be much higher than 2po). 
The second factor is always positive. Consequently the sign 
of AE sym (p, T) is determined by the third factor. If R m /R T >1 
the symmetry energy decreases with temperature, otherwise it 
increases with temperature. The splitting of the neutron and 
proton EFM influences the AE sym (p, T) explicitly through R m 
and indirectly through R T . Now let us investigate how the R T is 
eventually also influenced by the EFM splitting. Before com- 
ing to this point, we show in Fig. [2] the systematic calculated 
results with Skz-1, Skz2, and Skz3, respectively. We notice 
that the predictions from Skz-series Skyrme interactions are the 
same for the isoscalar part but different for the isovector part of 
the equation of state and moreover the Skz-1, Skz2, and Skz3 
belong to groups IV, V, VI, respectively. Therefore they are 
very suitable for exploring the correlation between the temper- 
ature dependence of the symmetry energy and the nucleon EFM 
splitting. In Fig. [2] the top panel shows the density dependence 
of symmetry energies at r=0, 5, 10, 20 MeV, respectively; the 
second panel shows the energies per nucleon in NM and SM at 
r=0 and 20 MeV, where the energies per nucleon in SM are the 
same for three interactions because they have the same isoscalar 
part; The third panel shows the nucleon(neutron) EFMes in SM 
and NM, and the bottom panel shows the density dependence 
of R,„, R T , and R,„/R T at T-20 MeV, respectively. The results 
calculated with other Skz-series interactions are similar to the 
one of these three Skz interactions. We first investigate the cal- 
culation results with Skz-1, the strong neutron and proton EFM 
splitting case with /?",=1.84, one sees from Fig.2(cl) that the 
neutron EFM in NM is larger than 1 and increases with density 



3 




0.5 1.0 0.0 0.5 1.0 0.0 0.5 1.0 1.5 



P'Pn 



P'Pn 



P'Pn 



Figure 2: (color online) The symmetry energy at T=0, 5,10,20 MeV, the energy 
per nucleon at T=0 and 20 MeV , the neutron and proton EFMes , the R m , R T , 
and R„,/R T at 7=20 MeV in NM and SM as function of density calculated with 
Skz- 1 , Skz2 and Skz3 interactions, respectively. 



quickly, the proton EFM is smaller than 1 and decreases with 
density even faster, and the nucleon EFM in SM is just between 
them and decreases with density, which leads the R m to increase 
with density quickly as shown in Fig.2(dl). From Fig.2(bl), 
one sees that the energy per nucleon in NM increases with den- 
sity and that in SM decreases with density for both T=20 and 
T=0 MeV cases but the increasing slopes in NM and decreasing 
slopes in SM with density are different for T-2Q and T=0 MeV. 
As a result, the symmetry energy at T=20 MeV becomes higher 
than that at T=0 MeV when p >0.75po(Fig.2(al). It is more ev- 
ident from Fig.2(dl). For T=20 MeV case, the R m /R T becomes 
smaller than one when the density exceeds about 0.75po and 
therefore the symmetry energy to be larger than that in T=0 
MeV case from Eq.(9). In the middle and right columns, we 
show the results calculated with Skz2 and Skz3. For Skz2( 
m* > to* with#°=1.13)andforSkz3(m* < to* with #£=0.91), 
the splitting of the neutron and proton EFM is small and both 
neutron and proton EFM are close to 1 . Accordingly, the R W /R T 
becomes larger than 1 at T=20 MeV for both Skz2 and Skz3 
cases and thus the symmetry energies for those two cases are 
all decreasing with temperature in whole densities. The extent 
of decreasing depends on the magnitude of R m /R T . The cal- 
culation results for v-series Skyrme interactions have similar 
behavior as those of Skz-series. 

To understand what factors will influence the R T , let us start 
from the expression (0 or (0 for n q (p). For cold matter, the 
n q (p) is a step-function, the momenta of all particles are within 
the Fermi momentum sphere with k„f = k p f = (3n 2 p/2) 1 ^ for 
SM and k„f = (3n 2 p) [/i for NM. For hot matter we know from 
expressions @ or (D that the deviation of n q (p) from a step- 
function depends on the density (related to the chemical poten- 
tial term in (0 and (0) and nucleon EFM in the matter(related 
to the kinetic energy of nucleon in (0) when temperature is 
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Figure 3: (color online) The neutron(proton) occupation number distribution 
functions in NM and SM at T=0, 20 MeV calculated with Skz-1 and Skz2, 
respectively. 



fixed. More explicitly,the deviation of n q (p) for hot matter from 
a step-function is larger for a system with lower density and 
larger nucleon EFM. In Fig. [3] the left panel shows the n q (p) 
at p — 0.3po and the right panel shows the n q (p) at p — 1 .5po 
at T=0, 20 MeV calculated with Skz-1 and Skz2, respectively. 
The solid lines in the figure denote n q (p) for neutrons(protons) 
in SM and dashed lines denote n„(p) for neutrons in NM at 
T=0. The excessive kinetic energy needed for conversion of all 
protons in SM into neutrons in NM comes from the excessive 
neutrons in n q (p) in NM as compared with that in SM, which is 
one of the origin of the symmetry energy for cold matter. For 
the matter at finite temperature, the calculated results of n q (p) 
are shown by the lines with symbols in the figure. The calcula- 
tion results are just as expected. From Fig. 3(a) and Fig. 3(c) for 
low density p = 0.3po case we find that the n q (p) at T=20 MeV 
for both 6 = and 6=1 deviates from step-function largely 
but there is no obvious difference between the results calcu- 
lated with Skz- 1 and Skz2 because of low matter density. For 
high density p = 1 .5po case, the results of n q (p) (right panel) in 
SM (6=0) at T=20 MeV calculated with Skz-1 and Skz2 are the 
same because of the same iso-scalar part of two interactions and 
both deviate from step-function not largely. While, the results 
of n q (p) in NM (6=1) calculated with Skz-1 and Skz2 become 
much different. For the case with Skz2(Fig.3(d)), the deviation 
of the nucleon n q (p) in NM from a step-function is almost the 
same or even weaker than that in SM and the ratio R T becomes 
smaller than l(see Fig.2(d2)), which leads the R m /R T at T=20 
MeV to be larger than 1 because R m > 1 and consequently the 
symmetry energy at T=20 MeV becomes smaller than that at 
T=0. While for the case with Skz-1, the nucleon n„(p) at T=20 
MeV in NM deviates from that at T=0 MeV largely (Fig. 3(b)) 
and its very long tail means that a large portion of neutrons 
inside the Fermi momentum sphere are excited to the outside 
of the Fermi momentum sphere and quite large part of them 
have very large momenta. And consequently the kinetic energy 
increases largely. Thus, the R T becomes much large at high 
density and even larger than the R m which is larger than one. 
Finally, R m /R T at high densities becomes smaller than 1 and the 
symmetry energy at high density at T=20 MeV is enhanced as 
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compared with cold matter. From above discussion, we find 
that the splitting of the neutron and proton EFM influences the 
R T through the difference between the n q {p) in SM and that in 
NM. This effect on R T can even be amplified as compared with 
that on R„, for some Skyrme interactions. Thus, ultimately, the 
splitting of the neutron and proton EFM influences the trend of 
the density dependence of symmetry energy in hot matter. For 
some Skyrme interactions like Skz-1 which has strong split- 
ting of the neutron and proton EFM and m* > ;«*, the onset 
of the TrTDSE at a certain density, usually at supernormal den- 
sity, leads the symmetry energy to be stiffer in hot matter than 
that in cold matter. For others, the splitting of the neutron and 
proton EFM leads the symmetry energy to be more soft in hot 
matter than that in cold matter. Generally, if the Skyrme in- 
teraction belongs to group IV and V with its R® n being much 
larger than 1, the TrTDSE will occur at certain density. The 
density for the onset of the TrTDSE depends on the magnitude 
of the splitting of the neutron and proton EFM in asymmetric 
matter. Concerning experimentally probing the TrTDSE or the 
temperature and density dependence of the symmetry energy, it 
requires to simultaneously measure the observables such as the 
transverse momentum spectrum of protons(pions,or other light 
charged particles) to determine the temperature [22] and those 
to be sensitive to symmetry energy at high densities such as the 
ratio between the yields of different charged pions, if /jt + , to 
extract the symmetry potential at high densities in heavy ion 
collisions with isospin asymmetric nuclei at energies of hun- 
dreds to around 1 GeV per nucleon. At present, highly con- 
tradictive results about the symmetry energy at high density 
extracted from experiments may have some relation with the 
unclear temperature dependence of the symmetry energy. 

In summary, we have investigated the temperature and den- 
sity dependence of the symmetry energy in hot nuclear matter 
by means of the approach of the thermodynamics with Skyrme 
energy density functional. We first classify the Skyrme inter- 
actions into 7 groups according to the range of the calculated 
neutron and proton EFM varying with density in the NM limit. 
Then we study the temperature dependence of symmetry energy 
for different Skyrme interactions. We find that for some Skyrme 
interactions the symmetry energy decreases with temperature at 
whole density region and it means that the symmetry energy in 
hot matter is softer than that in cold matter. But for Skyrme in- 
teractions belonging to groups IV and V (for these two groups 
m* > m*), the TrTDSE phenomenon may occur at certain den- 
sity, which makes the symmetry energy at high density to be 
stiffer. In order to explore the correlation between the neutron 
and proton EFM splitting and the temperature dependence of 
the symmetry energy in hot nuclear matter we introduce the ra- 
tio R,„ between the nucleon EFM in NM and that in SM, which 
is directly related with the neutron and proton EFM splitting 
in asymmetric matter and the ratio R T between the kinetic en- 
ergy density increasing from T-Q to T = T in NM and that in 
SM. The effect of the effective mass splitting on R T is indirect 
but for some Skyrme interactions it can be amplified as com- 
pared with that on R m . Only for those Skyrme interactions with 
R m /R T <1 in group IV and V the TrTDSE phenomenon can oc- 
cur at high densities and it ultimately makes the trend of the 



density dependence of symmetry energy to be more stiff in hot 
matter. The onset density of TrTDSE depends on the magnitude 
of R® r The larger R®, is, the lower onset density is. Finally, we 
would stress that the symmetry energy in hot nuclear matter not 
only depends on symmetry potential part but also on the split- 
ting of the neutron and proton EFM in kinetic part to a certain 
extent. Therefore, the study of the temperature dependence of 
the symmetry energy in addition to the density dependence is a 
matter of significance for extracting the symmetry potential and 
also the splitting of the neutron and proton EFM in asymmetry 
nuclear matter. 
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